Condensate density and superfluid mass density of a dilute Bose gas near the 

condensation transition 
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We derive, through analysis of the structure of diagrammatic perturbation theory, the scaling 
behavior of the condensate and superfluid mass density of a dilute Bose gas just below the con- 
densation transition. Sufficiently below the critical temperature, T c , the system is governed by the 
mean field (Bogoliubov) description of the particle excitations. Close to T c , however, mean field 
breaks down and the system undergoes a second order phase transition, rather than the first order 
transition predicted in Bogoliubov theory. Both condensation and superfluidity occur at the same 
critical temperature, T c and have similar scaling functions below T c , but different finite size scaling 
at T c to leading order in the system size. Through a simple self-consistent two loop calculation we 
derive the critical exponent for the condensate fraction, 2/3 ~ 0.66. 
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The behavior of a dilute interacting homogeneous Bose 
gas at the condensation temperature, T c , and in its neigh- 
borhood, is a delicate critical problem [1-3, and references 
therein]. Even though in a dilute gas, where na 3 <C 1, 
with n the number density, and a the s-wave particle- 
particle scattering length, perturbation theory fails close 
to T c , and fluctuations cannot be neglected; e.g., al- 
though the shift in T c from the free gas transition tem- 
perature, T°, is linear in a, the first order perturbation 
result vanishes, while all higher order terms are divergent. 

We calculate here the dependence on a of the conden- 
sate density, n , and the superfluid mass density, p s , of a 
dilute Bose gas, for temperatures T just below T c , where 
t = (T c - T)/T c is of order a/A. Here, A = {2n/mT) 1 ' 2 
is the thermal wavelength in units with K = 1, and m is 
the particle mass. As long established, the Bogoliubov 
(mean field) approximation fails close to T c ; it leads to 
a first order phase transition (e.g. 0]), vs. the second 
order transition expected for the universality class of the 
Bose gas. A similar phenomenon occurs in relativistic 
(j) 4 theory pj. We derive the general scaling structure of 
both no and p s in the critical region, as functions of t 
and a j A, which connects the mean field solution to the 
critical behavior of a continuous phase transition. The 
scaling functions for no and p s are similar, and imply 
that in the (dilute) interacting Bose gas the phenomena 
of condensation and superfluidity occur at precisely the 
same temperature. A further consequence is that in the 
very dilute limit, a — > 0, where the shift of the critical 
temperature is linear, T c — Xj? ~ a, {ppl, the conden- 
sate and superfluid mass densities at the ideal gas critical 
temperature, Tj?, both vary linearly with a. We also cal- 



culate the leading order finite-size corrections to no and 
p s at T c , which provides insight into the difference of the 
numerical results of Refs || and Q for the shift of the 
critical temperature. 

Our approach is to study the particle densities via di- 
agrammatic perturbation theory, working to the lowest 
needed order in a and no. To deduce the behavior of 
the superfluid mass density we employ Josephson's rela- 
tion between p s , no, and the long wavelength limit of the 
single particle Green's function [0J|J^1. 

The particle density, n, in the condensed phase is a 
function of a, no, and T, and has the form, n(a, no, T) — 
no + n(a, no, T), where h(a, no, T) is the density of non- 
condensed particles (with momentum k ^ 0). At the 
transition temperature, n(a, 0, T c ) — n; thus writing 
An = n(a,no,T) — n(a, 0,T), we have no + An = 
n(a, 0,T C ) — n(a, 0,T). Since the difference of T c from 
the ideal gas transition temperature, T°, is of order a pj, 
we may, to lowest order in a, for t of order a/ A, replace 
the difference on the right side by n(0, 0, T c )—h(0, 0, T) — 
n ((T c /T) 3 / 2 - 1)) ~ |ni. Up to corrections of order a 2 , 
we have then 
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Equation (Q) implicitly determines the condensate frac- 
tion, no(a,i), as a function of a and t in the critical re- 
gion, t Si a/A or no/n Ss a/A. 

It is simplest to calculate n(a,no,T) in terms 
of the matrix Green's function, G(rt,r't') = 
-i ((T(*(rt)*t( r 'f'))) _ (#t( r 't'))(#(rt))), where the 
two component field operator is ~^>{rt) = N)(rt),'ip'(rty\. 
The Fourier components of Q^ 1 have the form, 
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where the z n = 2-nnT are Matsubara frequencies (n = 
0, ±1,±2...), £ft = k 2 /2m, and the T,ij(k,z n ) are the 
corresponding self-energies. The chemical potential, fi, 
depends here on no- 

The non-condensate density, h, is then found from 



i(o, no, T) = -T 2 y ^Gn(* 
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with 
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Quite generally, An, to leading order in a and no, is given 
by the z n = contribution only: 



An = -T 



d 3 k 
(2^)3 



[G n (fc,0)-G(fc,0)], 



(5) 



where G(fc,0) = lim„ ^o Gn(fc, 0); integrated over /e, 
G(fc,0) determines the leading order shift of the critical 
temperature due to interactions [jjj. The Hugenholtz- 
Pines relation , 



/i = Sn(0,0)-Si 2 (0,0), 



specifies /1 as a function of no. In the zero fre- 
quency sector, Sn(fc,0) = £22(^,0) and £i 2 (fc, 0) = 
E 3 i(*,0), so that lim fc ^ [(/i - £n(fe,0)) (//- E a2 (A,0)) 
— Ei 2 (fc, 0)Ei2(fc, 0)] = 0, and thus the excitation spec- 
trum is gapless. In the following we drop the Matsubara 
frequency index, always referring to the zero frequency 
components. 

The lowest order mean field self-energies, En = £"'/ = 
2.9(^o + n )j E12 = E^ = gno, and u = g(no + 2n), where 
g = Ana/m, lead to the usual gapless Bogoliubov excita- 
tion spectrum. The mean field contribution to An, from 
Eq. §, is 
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Since the contribution of this term in Eq. (Q) is cx — n^ 2 
we find two possible solutions of (Q) for n at the mean 
field critical temperature, t m t = 0, namely n = and 
n = 47ra/A 4 ; intermediate values are not possible for 
t > 0. Thus Bogoliubov theory predicts a first order 
phase transition with a jump of the condensate density 
from to no = 47ra/A 4 0). However, as we discuss below, 



mean field can be valid only outside the critical region 
(where a/A <C |t| <C 1, and thus from Eq. (0), a <C n A 4 ), 
where it implies that no cx 

To go beyond mean field we analyze the structure of 
the self-energies by expanding En — Y,™/ and £12 — E" 2 ^ 
in a series in a and the mean field Green's functions G™ 
and Cf^/ , given by Eq. (j|) with the replaced by Ejy . 
We eliminate m in Q in favor of no using Eq. (^). How- 
ever, rather than using the gapless spectrum directly in 
a perturbative expansion, we write the propagators in 
terms of the mean field correlation length, £, as in [EJ, 
defined by 

M - (E^ - E™/) = (E n (0) - E 12 (0)) - (Ef/ - E™ 7 ) 



-l/2mC 2 
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Since the propagators remain formally infrared con- 
vergent we can derive the scaling structure of the self- 
energies by power-counting arguments. As above the 
transition, the ultraviolet part, when we neglect non-zero 
Matsubara contributions, has only a harmless logarith- 
mic divergence which can be removed by renormaliza- 
tion Q . The expansion of the self-energies beyond mean 
field starts at order a 2 ; furthermore E12 is formally at 
least of order no. Diagrams of order a K with k > 3 in 
the formal expansion contain vertices with two Green's 
functions entering; similar to the structure at T c , they in- 
volve the dimensionlcss combinations aC/A 2 and noA 2 £. 
The latter part originates from the dependence of Q m f 
(6) on 2mE^ ~ ano- Any diagram with an explicit power, 



p, of no can be generated from a corresponding diagram 
of power p — 1 in which a line is replaced by \/n at each 
of its ends. Thus each power of no involves one fewer 
three-momentum loop to be integrated over, replacing a 
structure of the form 2mT Jd 3 k/(k 2 + (~ 2 ) ~ C _1 A~ 2 , 
in a loop integral of order (a£/A 2 ) 2 , by no. The ex- 
plicit no dependences therefore enter in the combination 
«/A 2 ) 2 (n CA 2 ) = a 2 n C 3 /A 2 . Then with all momenta 
k scaled by we find the following scaling structure 
for the self-energies, 
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where the cry are dimensionlcss functions of dimension- 
less variables. In particular, for vanishing k, 

(E 11 (0)-E 12 (0))-(E™/(0)-E™ / (0)) 

= T^.s(g,n A 2 c), (10) 



where s is a dimensionless function. 
Equations (ft) and ([HI) imply that 
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C = — h(n Q X 4 /a), 
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where h is a dimensionless function. Then using Eqs. (|^) 
and (Jl|) in (§), we see that An has the scaling struc- 
ture An, = (a / X 4 ) f (noX 4 / a) . It immediately follows that 
close to T c , where no ~ a/ A 4 , the dimensionless func- 
tion / cannot be determined by a perturbation expan- 
sion in no A 3 or a/ A; therefore the predictions of mean 
field theory fail in this region. Finally from Eq. (|l|), us- 
ing n ~ 1/A 3 to lowest order, we derive the basic scaling 
result in the critical region, 
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where / is a dimensionless function. In the mean field 
limit, x — > oo, we must have f(x — > oo) ~ x, whereas the 
theory of critical phenomena implies a power-law behav- 
ior in the opposite limit, f{x — > 0) ~ x 2/3 , where /3 is the 
critical index for the order parameter, {ip) = ^/n Q . 

We see that for constant tX/a, no varies linearly with 
a. As a — > 0, T c varies linearly in a Q]; thus at the ideal 
gas transition temperature t = to = (T c — T c °)/T c ~ a, 
for a — ► 0, the condensate fraction varies linearly with a. 
This result is consistent with Leggett's weak variational 
bound that at to , no is bounded above by terms of order 

We now calculate the scaling function f(x) explic- 
itly within a simple model beyond the Bogoliubov ap- 
proximation. Introducing U{k) = 2m(En(fc) — En(0)), 
we have e k + En(jfe) - \i = (k 2 + U(k) + 2mY, 12 )/2m. 
Taking for £12 (fc) only the first order diagram in no, 



trnf 
J 12 



gno, we may write Eq. (pf) as 
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dk 



(k 2 + U{k)) {k 2 + U(k) + Amgn a ) ' 



(13) 



Were we to neglect U(k), we would derive the mean 
field result (Q); rather, we determine U(k) from a self- 
consistent two- loop calculation, as in [uj, 



U(k) 



-img 2 T 2 



d 3 p 



(2tt) 3 J (2tt) 3 
G(p + q)G(p) [G(k 



q)-G(q)}, (14) 



where G~ 1 (k) = /1 — Sk — U(k)/2m is the inverse of the 
(z n = 0) Green's function at the transition temperature. 
A free particle spectrum in G would leads to a logarithmic 
infrared divergence. As in the calculation of the critical 
temperature, self-consistency of U(k) at this level implies 
that it has the approximate structure, U{k) = k^ 2 k 3 / 2 , 
for non-zero k < k c , with k c = 32(27r/15) 1/2 a/A 2 w 
20.7a/A 2 §. A change of the power-law behavior of the 
free propagator to G[k) ~ —k~( 2 ~ 71 \ for k — > 0, with 
n > 0, is only possible precisely at T c , where the correla- 
tion length diverges. As long as no 7^ 0, the Bogoliubov 
operator inequality M, — Gn(fc,0) > mn /nfc 2 , does not 



permit n > 0. Therefore the z n = spectrum remains 
quadratic, n = 0, as k — > 0, everywhere below and above 
the critical temperature. 

However, for n A 4 a, the details of the spectrum 
at small k do not enter, and we may approximate the 
self-energy as: 



U(k) 



kl /2 k 3 ' 2 
k 2 



k < k c 
k > k c 



(15) 



which in Eq. (J131) leads to 



J (32n a/3fc c A 2 )ln(167rn a/fc 2 ) : n A 4 < a 
An ~ j -(2tt 1 / 2 /A 3 ) (n A 2 a) 1/2 : n A 4 » a. (16) 

The second line is the mean field result (g). Taken 
literally, this model would again predict a first order 
phase transition; however the logarithmic term indicates 
a change in the power law behavior close to the critical 
point. To determine this relation we invert Eq. (|12|), us- 



ing n ~ A~ 3 to note that in/no must be a dimensionless 
function of the variable A 4 no/a. In the limit no A 4 -C a, 
we may write, using the upper result in Eq. ([l6|), 

3 A 32a , / 167rnoaAA , „, 

no + An= 2 nt ^ no [ 1+ 3k-X^[-l^))^ 

which is the first term in an expansion of the scaling form 



nt ~ no 



(18) 



in the formal limit a/k c X 2 — > 0, consistent with our ap- 
proximation of U(k). Inverting, we derive n ~ t 2 @ , with 
2/3 = 1/(1 + (5/671-) 1 / 2 ) ~ 0.66, in excellent agreement 
with the value 2/3 ~ 2/3 expected for this universality 
class. In the other limit, noA 4 ^S> a this model calcula- 
tion simply approaches the mean field result, no ~ t. 

Below T c the condensed system is superfluid, with a 
superfluid mass density, p s , related to no by Josephson's 
sum rule @,|§, 
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k^o k 2 G n (k,0)' 
Using the explicit form (||) for Gn(k,0), we have 



d 

p s = nam + 2nom 2 ^-(Sii(fc) - T, 12 (k)) 



k=t) 



(19) 



(20) 



for T < T c . This result implies that precisely at T c , 
the superfluid fraction vanishes with the condensate frac- 
tion. Above T c , the superfluid density vanishes, as can 
be directly derived by calculating the transverse current- 
current correlation function p"^ . 

Further, we immediately see from Eqs. (|2^) and (|^) 
that the superfluid fraction in the neighborhood of the 



3 



ideal gas transition has the same scaling behavior as we 
found for the condensate density: 
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however the scaling function f p is in general different 
from the scaling function /. Reference (ll| obtained 
the scaling function for p s in the dilute limit to order 
e = 4 — d, where d = 3 is the spatial dimension. In 
the mean field limit, t ^S> a/ A, the lowest order self- 
energies are independent of k 2 , so that from Eq. j20| ) 
the superffuid mass density coincides with uq to order a. 
Thus f p (x — * oo) ~ x. In the critical region, however, 
f p (x — > 0) ~ x 7 , where 7 is the critical index for the su- 
pcrfiuid mass density. Josephson's scaling relation gives 
7 = 2/3 — 77^ ~ 2/3 for the critical index of the superfluid 
mass density, where v ~ 2/3 is the critical exponent of 
the correlation length Qj. Since our model calculation 
above does not include the correct k — > limit, to which 
p s (but not no) is sensitive, it is therefore not suitable 
for calculating p s reliably. 

Let us turn to understanding the behavior of no and 
p s in large but finite systems, of linear scale L. The con- 
densate density, tiq — n — h L at T c is non-zero for finite 
L, and is found to leading order from 



n = n n 



T 



d A k 

„ W 
d 3 k 



2ir/L 



G(k,0 



■G(k,C), 



(22) 



where G is the infinite size Green's function. Since at T c , 
G(k -> 0) = -2mCC n /k 2 -v, where C is constant, and 
ng° = 0, we find, 



AC 



(l + r/)A 2 i \ L J 



(23) 



to leading order, neglecting a numerical factor dependent 
on the particular geometry of the finite system. Joseph- 
son's relation should still hold inside the critical region 
of finite size systems, with the limit of zero wavevector 
replaced by k — ► 2it/L. With this relation we have 



^1/3 _^£_ = 
mn (1 



r])X 2 n 2 / 3 T c ° (l + ?y)(C(3/2) 2 / 3: 



(24) 



where the total particle number is given by N = nL . 
Note that both Eqs. (23) and ( pi| ) are valid independent 
of the diluteness of the gas. Equation (H) agrees well 
with the numerical values of Ref. [O] . Since the limit as 
a — ► of rj is non-zero for an interacting Bose gas, the 
formal a — > limit of Eq. ( pi] ) does not, however, agree 
with the ideal gas value, given by the same formula but 
with rj = 0. Therefore, the procedure of Ref. H, to ex- 
pand the finite-size scaling results directly around the 



ideal gas limit, is not completely justified; that there was 
a difficulty in this method was already suggested by the 
disagreement of the calculated value of T c — Tj? with later 
lattice calculations H which do not rely on this assump- 
tion. Nevertheless, use of finite-size scaling raises new 
strategies for explicit calculations Jl5[ . 
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